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Abstract Sequences of line-of-sight (LOS) magnetogranis recorded by the Michelson- 
Doppler Imager arc used to quantitatively characterize photospheric magnetic 
structure and evolution in three active regions that rotated across the Sun's 
disk during the Whole Heliosphere Interval (WHI), in an attempt to relate the 
photospheric magnetic properties of these active regions to flares and coronal 
mass ejections (CMEs). Several approaches are used in our analysis, on scales 
ranging from whole active regions, to magnetic features, to supergranular scales, 
and, finally, to individual pixels. We calculated several parameterizations of 
magnetic structure and evolution that have previously been associated with flare 
and CME activity, including total unsigned magnetic flux, magnetic flux near 
polarity inversion lines, amount of cancelled flux, the "proxy Poynting flux." and 
hclicity flux. To catalog flare events, we used flare lists derived from both GOES 
and RHESSI observations. By most such measures, AR 10988 should have been 
the most flare- and CME-productivc active region, and AR 10989 the least. 
Observations, however, were not consistent with this expectation: ARs 10988 
and 10989 produced similar numbers of flares, and AR 10989 also produced a 
few CMEs. These results highlight present limitations of statistics-based flare 
and CME forecasting tools that rely upon line-of-sight photospheric magnetic 
data alone. 
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1. Characterizing Photospheric Magnetic Evolution 

Evolution of magnetic fields in solar active regions (ARs) affects the heliosphere 
in many ways. In the low corona, sudden magnetic evolution in flares and coronal 
mass ejections (CMEs) — which arise on timescales of minutes to hours — can 
launch powerful disturbances into interplanetary space. More gradual magnetic 
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evolution, such as the effective diffusion and transport of active region magnetic 
flux over the photosphere, also drives evolution in heliospheric structure (affect- 
ing, e.g., the structure of the streamer belt and the heliospheric current sheet), on 
typical timescales of weeks and months. Consequently, studying the evolution 
of active-region magnetic fields is crucial to understanding the structure and 
dynamics of the heliosphere, beyond understanding solar activity itself. 

Quantifying evolution of magnetic fields in photospheric magnetograms, in 
particular, can provide insights into heliospheric evolution. Photospheric mag- 
netograms reveal cross sections of the large-scale structure of active regions, 
which extend from the solar interior out into the corona. The coupling of the 
photospheric field to the coronal field implies that magnetic evolution at the 
photosphere will induce evolution in the coronal field. Note, however, that the 
coronal magnetic field can evolve independently of the photospheric field; for 
instance, MHD instabilities that arise in the coronal field might trigger fiares or 
CMEs {e.g. Forbes, 2000). 

Here, we present analyses of several aspects of photospheric magnetic evo- 
lution in the three active regions (NOAA ARs 10987, 10988, and 10989) that 
crossed the solar disk during the Whole Heliosphere Interval (WHI) , as observed 
in line-of-sight (LOS) magnetograms recorded with the Michelson-Doppler Im- 
ager (MDI) instrument (Scherrer et ai, 1995) on the SOHO satellite. The data 
are described in Section 2, our analysis methods and results are discussed in 
Section 3, and we conclude by discussing the implications of these results in 
Section 4. 

2. Data 

2.1. Magnetograms 

Our analysis of photospheric magnetic evolution in the WHI ARs began with se- 
lection of magnetograms from the database of full-disk, line-of-sight, 96-minute- 
cadence. Level 1.8.2 magnetograms from MDI (http:/ /soi. stanford.edu/magnetic/ 
Levi. 8/). These magnetograms are formed either from a single measurement, 
nominally of 30 seconds, or by summing measurements over five minutes. The 
30-second magnetograms are much more common in the dataset than the five- 
minute magnetograms; and the former arc noisier. To estimate the noise levels 
in each, we fit Gaussians to histograms of field strengths (which are actually 
pixel-averaged flux densities) in one sample magnetogram of each type, and 
took the noise level as lata. Each sample magnetogram was smoothed with a 
2D, three-pixel boxcar prior to histogramming into 3 G bins. While we expect 
noise to dominate the weakest-field pixels, and would therefore prefer to fit 
only the "core" of each histogram (Hagenaar et ai, 1999), in practice using 
ranges smaller than ±30G resulted in poor fits. Results fitting this range of field 
strengths differed little from those derived by fitting the whole histogram. Whole- 
histogram (Tfit values were ~ 14G and ~ 8G for single and five-measurement 
magnetograms, respectively; the sample histograms used are plotted in Figure 
1. Given the scarcity of the five-minute magnetograms, we assume a single noise 
level of 14G for simplicity. 
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Figure 1. Histograms of pixel-average field strengths in single- (thick) and five- measurement 
(thin) magnctograms. Gaussian fits to the histograms can be used to estimate the noise levels 
in each, which arc ~ f4G and ~ 8G, respectively. 

The MDI instrument measures the LOS field strengtli [Blos] averaged over 
each pixel. Tracking the evolution of the radial magnetic field [Bn] requires esti- 
mating since only -Blos was observed. We therefore assumed the magnetic 
field was radial, and applied cosine corrections to the LOS field in each pixel, 
Bj{ = _Blos/ cos(7), where 7 is the heliocentric angle from disk center to each 
pixel. To compensate for foreshortening, triangulation was used to interpolate 
the En data — regularly gridded in the plane of the sky, but irregularly gridded 
in spherical coordinates (0, (f) on the solar surface — onto points (a;, y) corre- 
sponding to a regularly gridded Mercator projection of the spherical surface, 
following Welsch et al. (2009). This projection is conformal (and so locally 
preserves shape), which is necessary to ensure displacements measured by our 
tracking methods were not biased in direction. This reprojection distorts length 
scales, such that apparent displacements are exaggerated by a factor of the secant 
of the apparent latitude, which we corrected after tracking. 

Since the MDI magnetograph only measures the LOS component of the pho- 
tospheric magnetic field, we chose to analyze only magnctograms in which the 
target AR was within about 45° of disk center. The first AR visible on the 
disk during WHI was AR 10987 on 24 March 2008, and we started our tracking 
analyses with the magnetogram recorded at 14:23 UT. Shortly thereafter, AR 
10988 rotated onto the disk, followed by AR 10989. By 2 April, AR 10989 had 
no sunspots, but we analyzed magnctograms through 4 April, when AR 10989 
was nearly 45° from disk center, ending at 14:23 UT. The time ranges over 
which each AR was tracked are listed in Table 1. The top-most magnetogram in 
Figure 2 shows all three WHI ARs on the disk. The magnctograms in each of the 
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Table 1. Intervals over whieh photospheric magnetie fields of 
WHI ARs were tracked on the central disk. 



NOAA AR 


Tracking Start 


Tracking End 


10987 


2008 Mar 24, 14:23UT 


2008 Mar 30, 04:47UT 


10988 


2008 Mar 26, 12:48UT 


2008 Apr 01, 06:23UT 


10989 


2008 Mar 28, 22:23UT 


2008 Apr 03, 14:23UT 



subsequent rows of the same Figure show each WHI AR roughly one day after 
we began tracking it (left) and one day before we stopped tracking it (right). 

2.2. Flare and CME Activity 

Our primary goal is to relate properties of magnetic structure and evolution at 
the photosphere to energy release in the corona, in the form of flares and CMEs. 

Webb et al. (2011) used several sources of data to compile a comprehensive 
list of CMEs over the WHI, and to determine their source regions. We used 
their results to determine CME productivity of the WHI ARs during the time 
intervals shown in Table 3. Only AR 10989 was CME-productive while tracked; 
it produced two during its tracking interval (and a further two in the hours just 
outside the interval over which it was tracked). 

While CMEs are generally responsible for the strongest heliospheric distur- 
bances (Gosling, 1993), flares can also affect the heliosphere. As detailed below, 
each WHI AR only produced relatively small flares during the interval over 
which it was tracked. But even weak flares can affect the heliosphere: for instance, 
small hard X-ray (HXR) bursts are often associated with interplanetary Type III 
radio bursts (e.g. Christe, Krucker, and Lin, 2008a), which arise from electrons 
escaping into the heliosphere (where they have been observed in situ). 

To characterize flare activity in the WHI ARs, we first consulted records of soft 
X-ray (SXR) flare emission observed with the GOES satellites in the flare catalog 
maintained by NOAA (ftp://ftp.ngdc.noaa.gov/STP/SOLAR\_DATA/SOLAR_FLARES/ 
FLARES_XRAY/), and found that most of the relatively weak flares during WHI 
were not attributed to any source region, meaning this catalog did not accurately 
characterize each AR's flare activity. Next, we consulted Sam Freeland's Lat- 
est Events Archive (LEA: http://www.lmsal.com/solarsoft/latest_events_archive. 
html) which differences EUV images to determine flare-source locations. The 
LEA flare list identifled the source ARs for many events that lacked source ARs 
in the GOES catalog. In the three cases where the NOAA and LEA source 
attributions disagreed, manual inspection of LEA difference images suggested 
the LEA attributions were probably correct. Some flares in the LEA list were 
incorrectly attributed to either: i) remnant ARs from previous rotations that 
were due to rotate back to the flares' positions when they occurred, or ii) ARs 
farther West than the true source AR. Both errors might arise from using active 
region locations from Solar Region Summaries (SRSs), jointly prepared once 
per day by NOAA's Space Weather Prediction Center (SWPC; http://www. 



SOLA: ms_rev2.tex; 14 January 2013; 15:51; p. 4 



WHI Photospheric Magnetic Activity 



2Dai9-C'J- ?P 



AR 10937 



4€C 



800 ^Qua 
AR 1MB7 




100 



1H ZDD 

AR losae 



»o no 
AR 1KBB 




100 12a 14a IBO 1BD 3DD 3ID UD 



a2D MO sfio ssd eun ssd a4D 
AR 109B9 



420 

400 

I; 

I 



r 




AM 
















580 








sto 






. J 






m 


340 






\ : 


320 






•■".iJ 



1»1S014OieO1US00£2)0£'4O 



Figure 2. Top: This full-disk magnetogram, cropped at ±30°, shows all three WHI ARs (from 
right to left, 10987, 10988, 10989). The cropped magnetograms in each of the subsequent rows 
show each AR roughly one day after we began tracking it (left) and one day before we stopped 
tracking it (right). Axes' labels are pixels from MDI full-disk magnetograms. 
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swpc.noaa.gov/ftpdir/forecasts/SRS/README) and the U.S. Air Force, that were 
issued many hours before flare time. While source latitudes and longitudes in 
the LEA list appear accurate for the events we studied, AR source attributions 
from both the NOAA and LEA event lists contain errors. 

A more robust approach to automatic identification of flares' source active 
regions uses data from the RHESSI satellite (Lin et ai, 2002), which is capable 
of localizing the source regions of hard X-rays (HXRs) emitted in flares (as well 
as imaging HXR emission in higli-fluence flares) . Recent improvements in the al- 
gorithms used to identify microflares (short bursts of HXR emission in RHESSI's 
6-12 keV energy band; see, e.g., Christe et al., 2008b), which are invariably asso- 
ciated with active regions, has enabled their inclusion in the most recent RHESSI 
flare catalog (http://hessi.ssl.berkeley.edu/hessidata/dbase/hessijlare_list.txt). To 
flnd the source AR for each flare in the RHESSI catalog, disk positions for 
the WHI ARs from the daily SRSs were used to develop linear models of for 
each AR's apparent latitude and longitude as a function of time. This approach 
enabled extrapolating the positions of ARs both before they were visible near 
the east limb and recognized as new ARs, and after they passed beyond the west 
limb. Flares were automatically associated by proximity to the nearest predicted 
AR location. Visual inspection of each automatic association confirmed that 
this method correctly assigned above-the-limb flares to ARs which, though not 
visible, were at or near the limb. Lightcurves from GOES were then searched 
to identify the transient X-ray enhancement associated with each RHESSI flare, 
and the enhancement's peak intensity was recorded. 

Flares and source AR attributions are listed in Table 2, and plotted in Figure 
3, which shows the position of each RHESSI flare with respect to the solar disk, 
and its AR association. We note that most of the flares for AR 10988 occurred on 
the eastern disk, and many of the flares for 10989 occured above the east limb. 
Consequently, the LOS magnetograms from MDI cannot tell us much about the 
photosphcric magnetic structure of these ARs during their most flare-productive 
intervals. 

In Table 3, we summarize flare activity for each WHI AR separately for two 
time intervals: first, over the whole time range from appearance over the east 
limb to its disappearance around the west limb; and second, during the limited 
time interval that the AR was tracked. Not counting RHESSI flares for which a 
peak in the GOES lightcurve could not be automatically isolated, we found that 
ARs 10987, 10988, and 10999 produced 30, 44, and 28 flares, respectively, as they 
crossed the disk, and 9, 11, and 5 flares, respectively, while they were tracked 
across the central disk. Abramenko (2005) used a weighted, time-average of each 
AR's flare X-ray fluxes observed by GOES as a "flare index" to quantify AR 
flare activity; consistent with the definition of GOES flare classes, a power of 10 
difference in weighting was used between flare classes. Given such a weighting, 
the M- and C-class flares produced by AR 10989 imply that its average flare flux 
greatly exceeds that of the other two ARs during either of the intervals in Table 
3. 
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Table 2. RHESSI HXR bursts in the 612 keV channel during the WHI interval from the 
RHESSI flare list. Only those events with valid positions are shown in this list. Positions are 
listed in hcliographic coordinates (in degrees) if on the disk or in heliocentric coordinates 
(in arcseconds) if above the limb. The associated GOES class is also shown for each event; 
some events occurred during GOES data gaps. 
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Table 2. Continued from previous page; see above for more information. 



# Date Start Stop Peak Coordinates GOES Class Source AR 
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Figure 3. Positions and source ARs for flares identified by RHESSI. 
Table 2. Continued from previous page; see above for more information. 



# 


Date 


Start 


Stop 


Peak 


Coordinates 


GOES Class 


Source AR 


97 


2008 Apr 03 


04:56 


04:57 


04:57 


S10W64 


A8.1 


10988 


98 


2008 Apr 03 


05:55 


06:02 


05:57 


S09W65 


B1.2 


10988 


99 


2008 Apr 03 


06:02 


06:10 


06:07 


S09W65 


A9.9 


10988 


100 


2008 Apr 03 


07:35 


07:37 


07:36 


S10W63 


A4.2 


10988 


101 


2008 Apr 03 


20:00 


20:11 


20:04 


S06W77 


no data 


10988 


102 


2008 Apr 03 


20:15 


20:20 


20:17 


S13W41 


no data 


10989 


103 


2008 Apr 03 


20:20 


20:28 


20:23 


S13W42 


no data 


10989 


104 


2008 Apr 03 


22:09 


22:14 


22:11 


S06W75 


AS. 7 


10988 


105 


2008 Apr 03 


22:22 


22:24 


22:22 


S06W75 


A5.8 


10988 


106 


2008 Apr 05 


03:13 


03:17 


03:13 


S13W59 


A5.1 


10989 


107 


2008 Apr 05 


05:33 


05:37 


05:33 


S13W61 


A6.4 


10989 


108 


2008 Apr 05 


11:02 


11:16 


11:05 


S12W64 


B1.3 


10989 



3. Analysis of Magnetic Structure and Evolution 

We quantitatively analyzed magnetic evolution in these three ARs in several 
ways, which we present roughly in order of decreasing length scale: starting from 
whole-active-region measures of magnetic structure and evolution, we consider 
evolution on progressively smaller spatial scales, down to individual MDI pixels. 

3.1. Large Scale Structure and Evolution 

Barnes and Leka (2008), Leka and Barnes (2007), and Welsch et al. (2009) have 
noted that total unsigned radial magnetic flux [$] in an AR is strongly correlated 
with its flare productivity. 
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Table 3. Observed Flares and CMEs from WHI ARs, from limb to limb (top three rows), and 
while tracked (bottom three rows). A, B, C, and M denote GOES flare classes. (RHESSI flares 
during GOES data gaps are not included in these tallies.) CMEs within a few hours before or 
after the start or end times shown are included in parentheses before or after the frequency during 
the interval, respectively. The M- and C-class flares produced by AR 10989 imply that its flare 
index, a measure of AR flare productivity (Abramenko, 2005), exceeds those of ARs 10987 and 
10988. 



AR 


Start Obs. 


End Obs. 


A 


B 


C 


M Tot. Flares 


CMEs 


10987 


Mar 23, 23;39UT 


Apr 05, 11:17UT 


11 


19 




30 


0(1) 


10988 


Mar 23, 18:51UT 


Apr 05, 11:17UT 


24 


20 




44 


(1)4 


10989 


Mar 24, 19:09UT 


Apr 05, 11:17UT 


5 


20 


2 


1 28 


8(1) 


10987 


Mar 24, 14:23UT 


Mar 30, 04:47UT 


2 


7 




9 





10988 


Mar 26, 12:48UT 


Apr 01, 06:23UT 


4 


7 




11 





10989 


Mar 28, 22;23UT 


Apr 03, 14:23UT 


2 


2 


1 


5 


(1)2(1) 



We calculated the total unsigned estimated radial flux, [$ = J dA|i3fl|], in 
each of the WHI ARs over the time interval each was tracked, compensating 
for distortion in pixel areas due to the Mercator projection. These unsigned 
fluxes are plotted in Figure 4. The average values of unsigned flux for ARs 
10987, 10988, and 10989 were 1.6, 2.0, and 1.2 xlO^^ respectively, over the 
interval when each was tracked. To minimize the effects of noise, only pixels 
with unsigned field strength greater than 40G were included in the sums; this 
is approximately three times our assumed noise level of 14G. Assuming the 
uncorrelated, uniform per-pixel noise level on flux of (Jq ~ 14G(Aa;)^, where Ax 
is the width of MDI full-disk pixels, a formal estimate of the error in total flux 
would be (Tsum = \/^>40gO'o; this is typically < 10^° Mx, which is not visible 
on these plots. Dips in the plots arise from lower noise levels in magnetograms 
averaged over five measurements, and give a better estimate of the noise level. 
Also, the magnetic field is highly correlated from frame to frame at 96 minute 
cadences (Welsch et ai, 2009), implying frame-to-frame variations in total flux 
also indicate noise levels. Note that MDI has a known problem with saturation 
in strong-field regions (Liu, Norton, and Scherrer, 2007), such that flux values 
reported here might be lower than the actual flux (Wang et ai, 2009). 

Based upon total unsigned flux alone, AR 10988 should have been the most 
active, and AR 10989 the least. 

Schrijver (2007) argued that strong magnetic fields concentrated near polarity 
inversion lines (PILs) of the photospheric LOS magnetic field are closely associ- 
ated with the occurrence of large flares. He developed a method to quantify the 
total unsigned magnetic flux [R] near strong-field polarity inversion lines (SPILs) 
of ARs, by generating a weighting map of proximity to SPILs and summing 
the product of the weighting map and unsigned LOS magnetic flux. Because 
strong fields along LOS PILs are correlated with strong gradients in LOS fields 
across PILs (Falconer, Moore, and Gary, 2003), the R parameter is quantitatively 
related to the length of strong-gradient PILs, which has been found by Falconer, 
Moore, and Gary (2003) and Falconer, Moore, and Gary (2006) to be associated 
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Figure 4. Total unsigned flux [^J versus time in the three WHI ARs, as each crossed the 
central solar disk. Dips arise from lower noise levels in magnetograms averaged over five 
measurements. Larger values of "I" have been associated with greater flare activity. 



with CMEs. Barnes and Leka (2008) and Welsch et al. (2009) also found R to 
be associated with flare activity. 

Emulating Schrijvcr's approach here, we computed R values for each of the 
WHI ARs as they crossed the disk, choosing (as he did) a strong-field threshold 
of 150G and FWHM of 15Mm in the Gaussian used in the convolution to 
compute the SPIL weighting map (see Welsch and Li, 2008 for more details 
of the procedure). Our results are plotted in Figure 5. Note that our values 
of R are in units of Mx; Schrijver's were in units of G, summed over the 
weighting map. The MDI disk-center pixel area of ~ 2.2 x 10^® cm^ is an 
approximate conversion factor. Note also that Schrijver (2007) apparently used 
MDI Level 1.8.0 magnetograms, while we use Level 1.8.2 magnetograms (see 
http://soi.stanford.edu/magnetic/Levl. 8/). Since absolute flux densities in the 
latter are higher than in the former by a spatially varying factor of approximately 
1.6, our R values cannot easily be compared with Schrijver's. As with estimating 
total unsigned flux, uncertainties in R computed assuming an uncorrected, 
uniform per-pixel noise level [uq] are are too small to be visible on the plot. 
Instead, we have crudely estimated uncertainties from the R values themselves 
by computing the standard error in the mean with a five-point boxcar window. 
This simplistic approach is reasonable because, as noted by Welsch et aZ.(2009), 
changes in magnetic fields are relatively minor over a few hours, so successive 
measurements of R can be interpreted as repeated, independent measurements 
of nearly the same physical value. 

As with total unsigned flux, expectations of flare activity based upon R values 
would suggest AR 10988 should be most active, and AR 10989 should be least 
active. 



SOLA: ms_rev2.tex; 14 January 2013; 15:51; p. 11 



B.T. Welsch, J.M. McTiernan, and S.D. Christe 




3.2. Intermediate-Scale Evolution: Feature Tracking 

Automated tracking of "features" in magnetogram sequences has been used to 
understand processes governing evolution of the photospheric magnetic field, 
including flux emergence and dispersal [e.g. DeForest et al., 2007 and references 
therein). Features in magnetograms have been identified in several ways, for 
instance by: "clumping" collections of contiguous, like-polarity, above-threshold 
pixels {e.g. ParncU et al., 2009); "downhill" segmentation, which identifies dis- 
tinct "hilltops" in absolute field strength [e.g. Welsch and Longcope, 2003); and 
curvature, identifying the convex cores of hilltops in absolute field strength {e.g. 
Hagcnaar et al., 1999). Tracking involves the association of identified features 
between successive frames. 

Feature tracking can be used to identify episodes of flux cancellation, in 
which closely spaced, opposite-polarity features simultaneously appear to lose 
flux (Livi. Wang, and Martin, 1985). Since flux cancellation has been observed 
in prominence formation (Martin, 1998), shown to increase magnetic free energy 
in some circumstances (Welsch, 2006) , and proposed as a CME initiation mech- 
anism {e.g. Linker et al., 2003), we tracked downhill features in magnetogram 
sequences of the WHI ARs, to investigate cancellation rates in each AR. 

Using the YAFTA feature-tracking code (Welsch and Longcope, 2003, De- 
Forest et al., 2007; software and documentation are online at http://solarmuri. 
ssl.berkeley.edu/~welsch/public/software/YAFTA/) wc only included pixels above 
35G (a 2.5(7 threshold) in features, and required each feature to have at least four 
pixels, and at least one 70G pixel (a 5a peak threshold). We also required features 
to either persist for at least four frames, or interact with (fragment from, or merge 
with) a feature that persisted for four frames. The outlines of features identified 
in two successive magnetograms from AR 10987 are plotted over grayscale images 
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of the magnetic field in Figure 6. The outlines are color coded by feature label, 
and features have been matched between these magnetograms. Different-color 
outlines imply distinct features were identified; same-color outlines generally 
imply matching features between frames, although the color palette used is 
limited, so spurious color matches arc possible. While substantial evolution has 
occurred over the 96 minutes between these magnetograms, 90% of features in 
the first magnetogram were identified in the second, which is a typical survival 
rate. 

Having tracked the three ARs' magnetogram sequences, we then attempted to 
quantify canceling features in each. We first found opposite-polarity features in 
close proximity by looping over positive features, applying IDL's dilate function, 
with a (3 X 3) structuring element, to each feature's pixels; and seeking instances 
of overlap between the dilated feature's pixels and pixels of negative features. 
This essentially searched for negative features in nearest-neighbor pixels on the 
periphery of each positive feature. Hence, our definition of close proximity for 
defining cancellation is MDTs pixel size, ~ 2" ~ 1.4 Mm. For a pair of close 
features to have partially cancelled over a time step, we further required the 
features' centers of flux to approach, and the fiux of both features to decrease 
— and by doing so, we have implicitly assumed that neither feature disappears 
completely over a cancellation step. 

This approach is subject to many sources of error. Thresholds imposed on 
field strengths in features imply that some magnetic flux is not included in any 
feature. Features' positions jitter, and their fluxes fiuctuatc, both from noise 
and true evolution. Further, along the polarity-inversion lines where cancel- 
ing features are typically identified, each feature often interacts with several 
neighboring like-polarity features, and possibly other opposite-polarity features, 
perhaps exchanging flux or canceling with those features, too. Hence, there are 
many essentially random sources of changes in features' fluxes and locations. 
With higher-cadence data, additional constraints might be imposed, e.g. meet- 
ing our cancellation criteria over several time steps. However, a more-detailed 
analysis of cancellation in higher-cadence data lies outside the purview of this 
study. Further, we are primarily interested in comparing cancelled flux between 
the WHI ARs, not establishing absolute amounts of flux cancelled. 

Recognizing that our approach to quantifying cancelled fiux is subject to 
large uncertainties, we computed estimates in three ways: i) first, for each pair 
of canceling features identified, we took the minimum flux loss from the pair as 
the cancelled fiux in that event, and summed these fluxes from all cancellation 
events; ii) we limited our analysis to pairs of features that were observed to 
cancel over at least two steps, and took the average flux losses from such events; 
and Hi) again restricting our attention to features that were observed to cancel 
over at least two steps, we summed the minimum fiux lost from each pair in 
each event. The third approach is presumably the most restrictive method of 
estimating cancelled flux. 

For each AR, our tabulation of the total number of cancellation events, the 
number of multi-step cancellation events, and estimates for cancelled flux are 
listed in Table 4. Variations between the estimates for a single AR provide a 
measure of the uncertainties, although as noted above, the method used for the 
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Figure 6. The outlines of features — "iiills" in absolute field strength — identified in two 
successive magnetograms from AR 10987 are plotted over grayscale images of the magnetic 
field (white is positive polarity, black negative; saturation is set at ±250G). The outlines are 
color coded by feature label, and features have been matched between these magnetograms. 
Different-color outlines imply distinct features were identified; same-color outlines generally 
imply matching features between frames. 
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Table 4. Estimates of flux cancelled in the WHI ARs, computed three ways (see 
text); cancelled flux values are in units of [10^" Mx]. 



AR 


# Events 


# Multi-step 


All events, 
J2 min(<5'I>) 


Multi-step, 
E avg(5*) 


Multi-step, 
E min(<5'I>) 


10987 


58 


10 


5.0 


4.0 


1.9 


10988 


110 


21 


8.0 


8.1 


3.4 


10989 


42 


5 


2.9 


1.6 


0.8 



third estimate is the most restrictive. If our identification of canceUng features 
were essentially random, and mean feature fluxes in all three ARs were the same, 
then our estimates of cancellation events and cancelled fluxes would follow the 
ratio of unsigned magnetic fluxes in the ARs, 1.4:1.7:1 (see Section 3.1 above). 
While the numbers in each column do follow the ordering of these ratios, there 
is an excess of both cancellation events and cancelled flux in AR 10988 beyond 
that expected from variations in AR flux alone. 

As with both total unsigned flux and unsigned flux [R] near SPILs, our 
analysis of flux cancellation in the WHI ARs suggests that AR 10988 should 
have been the most active, and AR 10989 the least. 

3.3. Pixel-Scale Evolution 

3.3.1. Estimating Velocities 

According to Faraday's law, evolution of the radial magnetic fleld at the photo- 
sphere is governed by the curl of the electric fleld there, 

dtBr = -c(V X E)r . (1) 

Assuming that the electric fleld is ideal (equivalently, that the plasma's conduc- 
tivity is inflnite) implies E = — (v x B)/c, so evolution of the radial magnetic 
fleld can be related to horizontal variations in the velocity and magnetic flelds 
by the ideal induction equation, 

dtBr = [V X (v X B)],. (2) 
= -V,, • [whBr - . (3) 

The radial flux Sr of magnetic energy across the photosphere (the Foynting flux) 
depends upon E, and in the ideal approximation v, 

Sr = c[E X B],/47r = [vrBl - (v,, • B,0Br]/4^, (4) 

as does the rate of change of relative magnetic helicity (Berger and Field, 1984) 
in the corona dH/At due to the helicity flux across the photosphere, 

ATI r 

— = 2 y dA [(AC • B,)z;. - (Af • v,)i3.], (5) 
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where: the integration runs over the photosphere in the region of non-zero v 
and B; is the vector potential of the current-free (and therefore curl-free, or 
"potential") magnetic field that matches the observed photospheric radial 
field Br, i.e. V/, x = ; and = = V;, • A^. 

Demoulin and Berger (2003) argued that the "footpoints" of magnetic fields 
anchored in the photosphere appear to move in magnetograms with an apparent 
footpoint velocity u, which is related to the plasma velocity v by 

U = Vh- {Vr/Br)Bh . (6) 

They further argued that tracking methods applied to magnetograms, such as 
local correlation tracking (e.g. Chae, 2001), would estimate u, not Vf^. Welsch 
(2006) referred to u as the fiux transport velocity; uB^ has units of a flux 
transport rate (maxwells per unit length per unit time). Note that Equation (6) 
is a matter of definition, and so is made without approximation. In terms of u. 
Equation (3) can be written in a form analogous to the continuity equation, 

dtBr + V ■ (uBr) ^0 . (7) 

The statement that tracking methods will accurately estimate u, however, 
is a testable assertion; Schuck (2008) argued, in fact, that some "optical flow" 
methods are relatively insensitive to radial flows, so will return a biased estimate 
u of the true u. Tests of several velocity estimation methods conducted by Welsch 
et al. (2007) using synthetic data from MHD simulations, in which the true 
velocities were known, demonstrated that such methods were far from perfect, 
but that their flow estimates were highly correlated with the true flows. 

In terms of the flux transport velocity, the Poynting and helicity fluxes in 
Equations (4) and (5) become 

Sr = -(B,, •u,,)B,/47r (8) 
dH 



dt 



-2 J dA{A^ ■Uh)Br . (9) 



Note that estimating the Poynting flux requires knowledge of the horizontal 
magnetic field. 

The flux transport velocity can be estimated from magnetogram sequences, 
with finite cadence At and pixel size Ax, by applying a finite difference approx- 
imation to Equation (3), 

AB^ A • (uBr) , , 

where A is the spatial finite-difference operator. If a typical flow speed is wq, 
then the pixel-crossing time is tq = v^Ax. If At < tq, then Equation (10) 
approximates Equation (7); if, however. At > tq, the approximation can fail. 
This is analogous to the Courant-Friedrichs-Lewy (CFL) condition required for 
numerical stability in computational fluid dynamics. 

Welsch et al. (2009) used Equation (10) as the basis for applying two tracking 
methods, Fourier local correlation tracking (FLCT: Welsch et al., 2004, Fisher 
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and Welsch, 2008) and the differential affinc velocity estimator (DAVE: Schuck, 
2006), to sequences of 96-minute MDI full-disk magnetograms (similar to the 
current data, but Level 1.8.0 instead of Level 1.8.2). Both methods employed 
similar windowing parameters (eight and nine pixels for FLCT and DAVE, 
respectively) to localize the magnetic data used in estimating the flow at a 
given pixel, which probably has the effect of averaging over smaller-scale flows. 
Correlation coefficients between the field-weighted FLCT and DAVE flows' x— 
and J/— components were > 0.7 over the set of all the flows computed by Welsch et 
aZ.(2009), implying consistency between flow estimates using distinct methods. 
Mean and median speeds for FLCT were ~ 65 m sec^^, and for DAVE were 
~ 105 m sec~^; evidently, DAVE speeds are higher than those estimated with 
FLCT. Given the ~ 1.4 Mm disk-center width of MDI full-disk pixels, these 
speeds imply pixel-crossing times of '-^ 4 — 6 hr. Welsch et at (2009) also found 
that flow directions remained significantly correlated over similar time scales. 
Both of these facts imply that 96-minute cadence, full-disk MDI data can be 
used with Equation (10) as a valid approximation of Equation (7). 

We applied FLCT to the 96-minute cadence MDI full-disk magnetograms 
within 23 Mar 2008, 00:00:OOUT - 4 Apr 2008, 00:00:OOUT from the WHI 
interval, with a windowing parameter of eight pixels, an absolute field-strength 
threshold for tracking of 50 G, and a low-pass spatial wavenumber roll-off param- 
eter of 0.25 (see Fisher and Welsch, 2008). Pixels out to 45° from disk center were 
tracked. This choice of windowing parameter corresponds to ~ lOMm. a scale 
much larger than granules, but slightly smaller than supergranules (Hagenaar, 
Schrijver, and Title, 1997). 

Data cubes covering each of the WHI ARs' magnetic fields and their corre- 
sponding flows were then extracted from the full-disk magnetograms and inner- 
disk flow maps for the time intervals listed in the bottom three rows of Table 
3. Figure 7 shows two successive FLCT flow maps from the central, strong-field 
portion of AR 10987, from the same magnetograms shown in Figure 6. Some 
evolution in the flow pattern can be seen over the 96 minutes between the flow 
maps; but flows in many regions appear qualitatively similar at both times. 

3.3.2. Proxy Poynting Flux 

Welsch et al. (2009) analyzed relationships of flare activity with photospheric 
magnetic field and flow properties in a sample of a few dozen ARs. Of several 
quantities that they investigated, they found Schrijver's R value and a quantity 
proportional to ^{uBj^), where u = |u|, and the sum runs over the AR mag- 
netograms, to be most strongly associated with flare activity. Because uB^ has 
the same dimensions as the Poynting flux of magnetic energy (erg cm~^ sec""'^), 
Welsch et al. (2009) called this quantity the proxy Poynting flux. For consistency 
with Equation (8), we deflne the proxy Poynting flux here to be Su = ("-S|f)/47r. 
(Note that this definition differs from that of Welsch et a/. (2009), who defined 
Sr to include a sum over pixels — a total flux, instead of a flux density — and 
did not include 47r.) Li et al. (2010) estimated the proxy Poynting flux in AR 
8038 around the time of the well-known 12 May 1997 flare/CME, and found 
the cumulative flux reached ~ 10^^ erg over the four-day interval before the 
eruption. 
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Figure 7. Two successive FLCT flow maps from tlie central, strong-field portion of AR 10987, 
from the same magnetograms shown in Figure 6. Background grayscale shows estimated radial 
field strength, white showing positive flux, black negative, with saturation set at ±250G. 
White/black contours correspond to ±100, ±300, and ±500G. For clarity, only every other 
flow vector is plotted. Some evolution in the flow pattern can be seen; but flows in many 
regions appear generally similar at both times. 
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What is the physical basis for the relationship between the proxy Poynting flux 
and flare activity? Leaving aside the likelihood that estimates of Br and u are 
imperfect, in principle, 5*^ corresponds to part of the horizontal Poynting flux 
[Sh] (additional contributions arise from terms containing B^). The magnetic 
energy that powers flares and CMEs enters the corona from the solar interior, 
as an outward, radial Poynting flux [Sr]- It is plausible that there should be a 
correlation between 5*^ and flare/CME activity, although the corona can store 
magnetic energy, so there might typically be a latency between the introduction 
of magnetic energy into the corona and its release in flares/CMEs. For instance, 
Longcope et al. (2005) and Schrijver et al. (2005) report latency times associated 
with active-region magnetic fields of ~ 24 hours. (And in quiet-Sun regions, 
filaments/ prominences, which are interpreted as current-carrying [non-potential] 
magnetic field structures, can persist in the corona for weeks.) It is also plausible 
that the radial Poynting fiux [Sr] is significantly correlated with [St] and/or the 
proxy fiux [Sr\ , but no quantitative theoretical or observational characterization 
of such a correlation has been presented. 

Using the FLCT fiows we estimated for each WHI AR, we computed the 
proxy Poynting flux in that AR over the time interval between each tracked 
magnctogram pair. We then multiplied by (reprojection- corrected) pixel area 
AA and the time interval At between magnetograms, and summed over all pixels 
to express our results in terms of ergs of energy transported. Results are plotted 
in the top panel of Figure 8; in the bottom panel, we plot the cumulative proxy 
Poynting fiuxes over time. As with R, we estimated uncertainties in the top 
panel from a running-boxcar computation of the standard error in the mean. 
Assuming a combined fractional uncertainty of 0.6 for the terms in the product 
uB^ in each pixel resulted in error bars smaller than those shown. Summing 
the uncertainties in the top panel in quadrature for the cumulative result shown 
in the bottom panel yields error bars about the size of the steps in the bottom 
panel. 

Unlike the true Poynting fiux, which can be positive or negative {i.e. radially 
outward or inward), the proxy Poynting fiux is only positive. To the extent that 
the proxy Poynting flux reflects the outward radial fiux of magnetic energy, the 
values in these plots suggest the order of magnitude of energy transport in the 
WHI ARs. 

As with other magnetic parameters of the WHI ARs, values of the proxy 
Poynting flux suggest AR 10988 should have been the most active, and AR 
10989 the least. 

3.3.3. Helicity Flux 

Magnetic helicity is approximately conserved even when magnetic reconnec- 
tion occurs, meaning that it is not readily dissipated in the corona (Berger, 
1984), and can accumulate there as twisted magnetic flux emerges from the 
interior (Pevtsov, Canfield, and McClymont, 1997). It is plausible that helicity 
accumulation in the corona plays a key role in triggering CMEs (Low, 2001). 

Velocities inferred from tracking can also be used to estimate the flux of 
magnetic helicity across the photosphere (e.g., Chae, 2001, Kusano et ai, 2002, 
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Proxy Poynting Flux in ARs 987-989 
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Figure 8. The proxy Poynting flux [Su] multiplied by pixel area [A A] and the time between 
magnetograms [At] (top) and its sum over time (bottom) are both plotted versus time for the 
three WHI ARs, as eaeh crossed the eentral solar disk. Error bars for the cumulative plot would 
be about as large as vertical step sizes. Large values of X]pix '^R have been associated with flare 
activity. The sum over Sji can provide an order-of-magnitude estimate of the magnetic-energy 
transport rates. 
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Pariat, Dcmoulin, and Berger , 2005). Using the FLCT flows we estimated for 
each WHf AR, we computed the hehcity flux in that AR over the time interval 
between each tracked magnetogram pair, separately using i) Equation (9) with 
a Green's function method to calculated A^, and ii) the method presented by 
Pariat, Dcmoulin, and Berger (2005). Means, medians, and linear fits between 
the two methods agreed to within 3% for AR 10987, 1% for AR 10988, and 
6% for AR 10989. In the top panel of Figure 9, we show the differential fluxes 
of helicity for the WHI ARs computed from Equation (9), with uncertainties 
again computed from a running-boxcar calculation of the standard error in flux 
estimates. (As above, formally propagated errors were too small to appear on 
these plots.) In the lower panel of Figure 9 we plot the cumulative helicity fluxes 
of each WHI AR; only every third error bar is plotted, for clarity. 

A characteristic magnitude for helicity is the flux in the system squared. In 
Figure 10, we plot the cumulative helicity fluxes of each WHI AR, normalized 
by the square of the mean magnetic flux in each AR over the time period it was 
observed. The unnormalized, cumulative, helicity fluxes for ARs 10988 and 10989 
are similar; but since the relative magnetic fluxes of the ARs 10987, 10988, and 
10989 are 1.4:1.7:1, normalizing by their squared fluxes implies the characteristic 
hehcity for AR 10989 is substantially higher than for 10988. 

We have reviewed several magnetic parameters of the WHI ARs that might 
be correlated with flare and CME activity. Characteristic hehcity is the only 
such parameter for which AR 10989 has a larger value than the other ARs. 
The characteristic hehcity fluxes that we found, however, are far below unity, 
implying that the hehcity fluxes did not inject large amounts of twist into the 
WHI ARs, relative to their magnetic size. 

The total helicity present in an AR is the combination of that present in the 
AR when it emerged and the helicity injected after emergence. By extrapolating 
non-linear force-free flelds (NLFFFs) from a vector magnetogram, the instan- 
taneous helicity content of an AR can be estimated. Petri, Canou, and Amari 
(2011) have done this for the WHI ARs with SOLIS (Henncy et al, 2009) 
data, and report characteristic hclicities (see their Table 1) of 0.0055,0.0041, 
and 0.0017 for ARs 10987, 10988, 10989, respectively. Note that all of their 
characteristic helicity values are on the same order as those we report, but 
that all of their values are positive; and that their value of the characteristic 
helicity in AR 10989 is less than those of the other WHI ARs, undermining the 
hypothesis that characteristic helicity might explain variations in flare and CME 
productivity. 

As an aside, we note that some of our estimates of helicity changes also 
differ those of Petri, Canou, and Amari (2011) (see the right column of their 
Figure 23). They show AR 10987 with initially positive helicity that tends to 
decrease with time, which only partially agrees with our results, which show an 
initiafly negative helicity flux followed by a larger, positive hehcity flux. Our 
estimates show negative helicity steadily being added to AR 10988 over time, 
but Petri, Canou, and Amari (2011) found increasingly positive helicity over 
one time interval. Finally, the estimates of helicity changes in AR 10989 by 
Petri, Canou, and Amari (2011) arc relatively weak, and of mixed sign; but our 
results show a steady flux of negative helicity. 
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Figure 9. Top: the differential ffuxes of helicity for the WHI ARs, with uncertainties. Bottom; 
The cumulative helicity fluxes of each WHI AR; only every third error bar is plotted, for clarity. 
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Figure 10. The cumulative characteristic helicity fiuxes of each WHI AR: the cumulative 
hclicity flux, normalized by the square of the mean magnetic flux in each AR over the time 
period it was observed. The unnormalized, cumulative, helicity fluxes for ARs 10988 and 10989 
are similar, but the characteristic helicity for AR 10989 is substantially higher than for 10988. 

3.4. Inductive Flow Estimation 

Flows can be estimated from the induction equation alone, by first applying a 
Helmholtz decomposition to the horizontal electric field, 

cEft = -Vft0 + V,, X xf , (11) 

and using this result in Equation (1), resulting in a Poisson's equation for 

dtBr = Vlx ■ (12) 

If the electric field is assumed to be ideal, then cEh = — u x B^f , and a "Poisson 
flow" can be derived from x, 

uB, = -VftX • (13) 

We solved the finite-difference approximation to Equation (12) for the WHI 
AR magnetogram sequences using Fast Fourier Transforms (FFTs), and then 
derived Poisson fiows. Mean and median Poisson flow speeds for the pair of 
magnetograms shown in Figure 7 were 91 and 63 m sec~^, respectively, slow 
enough to be consistent with the finite-difference approximation. 

One problem with this approach is that fiuctuations in individual pixels due 
purely to noise or observational artifacts are ascribed to actual electric fields or 
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flows. In contrast, tracking methods such as FLCT and DAVE estimate flows at 
a given pixel using data from the neighborhood around that pixel (defined with 
windowing functions), and are therefore less susceptible to spurious fluctuations 
in field strengths in any individual pixel. 

Attempts to ameliorate this shortcoming by smoothing the field difference 
[ABji] rapidly degraded the consistency of the Poisson solution with the actual 
magnetic evolution. This consistency is checked using Equation (7) to calculate 
the expected magnetic evolution assuming uBji is known (e.g. Welsch et ai, 
2007). Temporal averaging of higher-cadence data, however, should counter- 
act the effects of spurious, pixel-scale magnetic fluctuations, but still enable 
consistency between derived flows and Equation (7). 

Unfortunately, statistical agreement between FLCT and Poisson-flow compo- 
nents is poor, with both linear and rank-order correlations below 0.1 for the 
magnetogram pair in Figure 7, both when weighted by and not. Proxy 
Poynting fluxes derived from tracking and Poisson flow estimates are signiflcantly 
correlated, but this could arise solely from spatial correlation of Bj^. Helicity 
fluxes derived from tracking and Poisson flow estimates are also poorly corre- 
lated. Unlike the proxy Poynting flux, helicity-flux density depends sensitively 
upon direction of flow, since u is dotted with A^. Given that FLCT and Poisson- 
flow components are weakly correlated, such disagreement in helicity fluxes is 
not surprising. 

Given these disagreements between results from each method, which flow esti- 
mate is more credible? Welsch et al. (2009) found flow estimates from the FLCT 
and DAVE methods — both of which arc optical flow techniques, but derive 
flow estimates in very different ways (Schuck, 2006; Fisher and Welsch, 2008) — 
were signiflcantly correlated. If weighted by Bn, linear correlation coefficients 
were > 0.7 for each flow component for more than 2000 magnetogram pairs 
in their dataset. In addition, Welsch et al. (2009) found that flow components 
derived from cither FLCT or DAVE persisted in time, with a decorrelation time 
of approximately six hours (see their Figures 9 and 10). In contrast, Poisson 
flows for AR 10987 have much lower frame-to-frame correlations than FLCT 
flows; they essentially decorrelate after one 96-minute interval. Both correlations 
between flows determined using the independent FLCT and DAVE methods and 
the persistence of FLCT and DAVE flows in time suggest that the FLCT results 
found here are more credible than the Poisson flows. 

Welsch et al. (2007) tested several velocity-estimation methods using synthetic 
magnetogram sequences extracted from MHD simulations, in which the flows 
were known. Later, Welsch and Fisher (2008) tested Poisson-flow estimates with 
the same synthetic data, and found them to be about as accurate as optical-flow 
estimates. Since the tracking and Poisson flows were both signiflcantly correlated 
with the known MHD flows, they were also significantly correlated with each 
other. The lack of agreement between FLCT and Poisson flows found here for 96- 
minute cadence, full-disk MDI magnetogram data suggests the Poisson method 
of estimating flows might fail when applied to such data. The Poisson technique 
might, however, fare better when applied to data with higher cadence and/or 
higher spatial resolution, such as data from either MDI in its high-resolution 
mode or HMI. Further study of Poisson flow estimates is planned. 
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4. Discussion 

By all of the statistical predictors of flare activity based upon magnetogram 
structure and evolution that we investigated, AR 10989 should have been signif- 
icantly less active than both AR 10987 and 10988. During the intervals that we 
tracked each AR, however, flare activity in AR 10989 was comparable to activity 
in the other ARs. Further, Webb et al. (2011) found AR 10989 to be much more 
CME productive than the other two ARs. 

Why did these statistical predictors fail in the case of the WHI ARs? Given 
that we only consider three ARs here, and that these predictors are purely 
statistical, we should not be too surprised that there are exceptions. Also, 
flares are rare events, and the time interval of our observations was limited, 
so drawing conclusions about activity levels from the small number of flares 
observed is problematic. While we only analyzed the ARs from the WHI for this 
study, stronger conclusions regarding relationships (or lack thereof) between 
ARs' photospheric magnetic properties and flares could presumably be made by 
studying a larger AR sample. Nonetheless, it would be useful if some physical 
basis for the exceptions we found could be identified. 

Unfortunately, the observations available to us are of limited utility in this 
regard. First, a practical concern: as noted above, we have only LOS mag- 
netogram sequences for the WHI ARs. Long-duration, regular-cadence vector 
magnetogram sequences (such as those available from HMI, or that might be 
available from a network of terrestrial, SOLIS-like instruments) can reveal much 
more about the structure and evolution of AR magnetic fields. 

Second, another practical concern: we are largely ignorant of the WHI active 
regions' magnetic structure and coronal activity over much of their lives. A 
significant component of AR 10989's activity occurred both before and after we 
tracked it across the central solar disk. Because the structure and evolution of its 
photospheric magnetic field could only be reliably estimated while on the central 
disk, we can only speculate about this field away from disk center. Welsch et al. 
(2009) showed that AR magnetic fields can evolve substantially over the course 
of 24 hours, and Longcope et al. (2005) and Schrijver et al. (2005) estimated 
coronal energy storage times on the order of a day. If the field of AR 10989 
evolved significantly between its time on the central disk and the pre- and post- 
tracking intervals, its photospheric magnetic structure during those more active 
invtervals might have been more consistent with statistical expectations based 
upon its coronal activity level. Similarly, we often only have hints about flares 
and CMEs produced by ARs beyond the limb. For instance, STEREO-A data 
show that on 05 April, AR 10987 at longitude W105 produced a CME just before 
16:00UT with a speed estimated from LASCO data at 962 km s^^, associated 
with an A-class GOES flare. (See Nitta's excellent compilation of data for this 
and other events at http://www.lmsal.com/nitta/movies/flares_euvi/) Burkepile 
et al. (2004) has shown that GOES flux is correlated with CME kinetic energy, 
with a linear correlation coefficient of 0.74. Similarly, the events analyzed by Qiu 
and Yurchyshyn (2005) show that GOES fiux is correlated with CME speed: we 
find linear and rank-order correlation coefficients of 0.79 and 0.54, respectively, 
in their sample of 13 events. These results suggest that the A-class flare observed 
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by GOES was most likely a partially occulted strong C-class or weak M-class 
flare. To address this ignorance, having more extensive magnetograph, X-ray, 
in-situ, and EUV coverage of the Sun would be useful. A satellite in co-rotation 
with the Sun could provide comprehensive observations of ARs from birth to 
death; and a constellation of satellites might provide full, "47r" coverage. 

Third, studying photospheric data alone neglects aspects of coronal magnetic 
structure that arc relevant to flaring. For instance, the large-scale structure of 
the coronal magnetic field, including perhaps the proximity of an AR to open 
magnetic fields, might play a key role in CME productivity. Sterling and Moore 
(2001) argued that open flux near NOAA AR 8210 enabled repeated, homologous 
ejections (see their Figure 3). Similarly, Liu (2007) reported that CMEs under- 
neath the streamer belt tended to be substantially slower than CMEs outside 
the belt, suggesting that the global magnetic environment around ARs can play 
a significant role in the CME process. Global-scale potential-field source-surface 
(PFSS) modeling by Petri, Canou, and Amari (2011) demonstrates that AR 
10989 lies almost completely outside the streamer belt (see their Fig. 21), perhaps 
easing the ejection of coronal fields into the heliosphere as CMEs. Differences 
in large-scale coronal magnetic environment might partially explain why fiare 
activity levels were broadly similar among the WHI ARs, but CME activity levels 
were not. (It is also true that fiares and CMEs appear to be distinct phenomena: 
very large fiares arc almost always associated with CMEs, but that the associa- 
tion is weaker with weaker fiares (Andrews, 2003).) It would be useful to study 
objective, quantitative measures of "proximity" between open fiux and ARs, 
and correlations between CME activity and such proximity measures. Hudson 
and Li (2010) noted that the flare-to-CME ratio increased this solar minimum 
relative to the last, implying the corona "was relatively easy to disrupt" during 
this minimum, also suggesting global properties of the coronal magnetic field 
might play some role in the CME process. 

Fourth, properties of photospheric magnetic fields might be essentially un- 
related to properties of the coronal magnetic field that actually generate flares 
and CMEs. Having analyzed a large sample of vector magnetograms, Leka and 
Barnes (2007) suggested "the state of the photospheric magnetic field at any 
given time has limited bearing on whether that region will be fiare productive." 
(Based upon the case studies here, we might add the evolution of the line-of- 
sight photospheric magnetic field over any given time interval also has limited 
bearing on whether that region will be fiare productive — but we have not 
analyzed enough cases to justify such a conclusion.) Efforts to measure properties 
of coronal magnetic fields using IR (SOLARC: Lin, Kuhn, and Coulter, 2004; 
CoMP: Tomczyk and Mcintosh, 2009; ATST: CargiU, 2009) and radio (FASR: 
Cargill, 2009) instruments could lead to breakthroughs in our understanding of 
fiare and CME physics. 

Finally, leaving aside all observational limitations, fiares might be unpre- 
dictable for more fundamental reasons, a possibility also noted by Schrijver 
(2009). 
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